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Recent quantum oscillation experiments in underdoped high temperature superconductors seem 
to imply two paradoxes. The first paradox concerns the apparent non-existence of the signature of 
the electron pockets in angle resolved photoemission spectroscopy (ARPES). The second paradox is 
a clear signature of a small electron pocket in quantum oscihation experiments, but no evidence as 
yet of the corresponding hole pockets of approximately double the frequency of the electron pocket. 
This hole pockets should be present if the Fermi surface reconstruction is due to a commensurate 
density wave, assuming that Luttinger sum rule relating the area of the pockets and the total number 
of charge carriers holds. Here we provide possible resolutions of these apparent paradoxes from the 
commensurate d-density wave theory. To address the first paradox we have computed the ARPES 
spectral function subject to correlated disorder, natural to a class of experiments relevant to the 
materials studied in quantum oscillations. The intensity of the spectral function is significantly re- 
duced for the electron pockets for an intermediate range of disorder correlation length, and typically 
less than half the hole pocket is visible, mimicking Fermi arcs. Next we show from an exact transfer 
matrix calculation of the Shubnikov-de Haas oscillation that the usual disorder affects the electron 
pocket more significantly than the hole pocket. However, when, in addition, the scattering from 
vortices in the mixed state is included, it wipes out the frequency corresponding to the hole pocket. 
Thus, if we are correct, it will be necessary to do measurements at higher magnetic fields and even 
higher quality samples to recover the hole pocket frequency. 



I. INTRODUCTION 

High temperature superconductors have been ad- 
dressed from a remarkable number of vantage points. 
Nonetheless many of the basic questions still remain un- 
resolved, and the notion of broken symmetries arising 
from a Fermi liquid has been traditionally discarded in 
favor of many exotic ideas. Here we revisit the Fermi 
liquid concept and a particular broken symmetry in re- 
sponse to a class of recent quantum oscillation exper- 
iments J-'^'^'i'^'^ We and othersii^i^iiaiiiii^iH have had 
some success in this respect. If these theories and ex- 
periments are correct, one will have to radically alter our 
twenty year-old view of these superconductors The to- 
tality of phenomenology cannot of course be explained 
without serious Fermi liquid corrections. But as long as 
the quasiparticlc residue is finite, wc hope that the low 
energy properties can be understood from our perspec- 
tive. 

Two paradoxes have arisen in the context of the quan- 
tum oscillation measurements. The first is the contrast 
between Fermi arcs observed in angle resolved photoe- 
mission (ARPES) experiments^^ on one hand and recent 
quantum oscillation experiments suggesting small Fermi 
pockets in underdoped YBaaCuaOg+s (YBCO)^i^i^i^'^ 
on the other. This is particularly clear in recent ARPES 
experiments where an effort was made to examine YBCO 
with similar doping as in the quantum oscillation mea- 
surements 

The second paradox is the non-existence of any evi- 
dence of the hole pockets in quantum oscillation mea- 



surements. Within the density wave scenario of wave 
vector Q = (tt, tt) (the lattice spacing set to unity), there 
should be two dominant frequencies in quantum oscilla- 
tions. One corresponding to the electron pocket at about 
500T and the other corresponding to the hole pocket at 
around 900 T. These are of course rough numbers cor- 
responding to approximately 10% doping, assuming that 
the Luttinger sum rule is satisfied in the mixed state, 
that is, the quantum oscillations reflect the normal state 
even if the measurements may lie within the mixed state. 

The concept of a broken symmetry is very powerful 
because deep inside a phase a physically correct effec- 
tive Hamiltonian can address many important questions, 
whereas our inability to reliably predict properties of 
even a single band Hubbard model, while widely pursued, 
has been a limiting factor. This is not an empty exer- 
cise, if new phenomena can be predicted or striking facts 
can be explained with some degree of simplicity. That 
broken symmetry both dictates and protects the nature 
of elementary excitations, determining the properties of 
matter, is important to emphasize. 

The suggested form of order, the d-density wave 
(DDW)^ explains numerous properties of these super- 
conductors, including the concomitant suppression of the 
superfluid densitj*i^. Hall numberfi^ and more recently 
the large enhancement of Nernst effect in the pseudogap 
state^^ in addition to the existence of a single-particle 
gap of dx2^y2 form above Tc- There are also theoretical 
reasons why DDW is a possibility. It competes favorably 
with other ordering tendencies in variational studies of 
extended Hubbard models with nearest-neighbor repul- 
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sion and pair-hopping termsi^i^ It is also realized in 
a class of two-leg ladder models with nearest-neighbor 
repulsion.— However, whatever form the correct Hamil- 
tonian takes, we know that it must favor d-wave super- 
conductivity (DSC). Such a Hamiltonian will almost cer- 
tainly favor DDW order as well, in light of the abundance 
of local Hamiltonians which do not discriminate between 
DSC and DDW order. In fact, two carefully designed dif- 
ficult polarized neutron scattering experiments have pro- 
vided tantalizing direct evidence of DDW order al- 
though other experiments have claimed otherwisej ^^'^^i^^ 
A natural enemy of the pristine properties of matter 
is disorder that is unavoidable in complex systems such 
as high temperature superconductors. The role of dis- 
order was emphasized in the original proposal of DDW 
order as a relevant competing order in the phase dia- 
gram of high-Tc superconductorSfii although our views 
of disorder have greatly evolved during the intervening 
years. It is this DDW order combined with disorder that 
would be the focus of the present manuscript in resolving 
the paradoxes stated above. The disorder considered are 
of two different types: (1) scattering due to impurities 
and defects and (2) scattering from vortices in the mixed 
states. 

We consider two kinds of intrinsic disoredr: (a) Gaus- 
sian white noise and (b) correlated disorder with a finite 
correlation length. In the momentum space the scatter- 
ing rate for correlated disorder will decay as exp(— g^Z^), 
where q is the momentum transfer between the initial and 
the final states, Id being the correlation length. There- 
fore, because of its smaller size, the states correspond- 
ing to the electron pockets arc scattered more than on 
the hole pockets. This is an interpretation of the phe- 
nomenon and is based on intrapocket scattering. An al- 
ternative interpretation involves the density of states on 
the Fermi surface of the electron pockets. In contrast, for 
white noise, scattering is independent of momentum and 
affects both pockets similarly. Disorder naturally has a 
strong effect on ARPES spectral function, which is sensi- 
tive to the coherence factors that are analogs of Wannier 
functions. In Shubnikov-de Haas oscillations it is only 
the averaged effect of disorder that enters by determin- 
ing the effective lifetime on the Fermi surface. Therefore, 
the role of disorder is quite different, as we shall explicitly 
see. 

Correlated disorder is also experimentally relevant. 
Unlike quantum oscillation experiments which probes 
bulk properties, ARPES is inherently a surface probe. 
In the relevant case of YBCO, as cleaved surfaces show 
that CuO and BaO terminations give different contribu- 
tions to the total photocmission intensity, with a hole 
doping Ufi = 30%, almost irrespective of the nominal 
bulk doping. This self-doping was controlled by evapo- 
rating potassium in situ on the cleaved surface, so as to 
reduce the hole content down to the value of underdoped 
bulk YBCO {S w 0.5) the doping level for which many 
quantum oscillation experiments are carried out. The 
potassium overlayer is likely to produce an effective cor- 



related disorder in the CuO plane. 

To explain the second paradox we shall adapt an anal- 
ysis of Stephen^S to include a normal state that exhibits 
DDW order. We shall see that the relativistic charac- 
ter of the nodal fermions of the hole pocket, as opposed 
to the nonrelativistic nature of the charge carriers of 
the electron pocket provides a possible explanation. If 
we denote the Dingle factor of the electron pocket by 
Ve = e^^^^"^'-', the Dingle factor of the hole pocket is 
Dfi w 2?g'^. This huge suppression may be the resolution 
of the missing frequency from the hole pocket. Here, 
LUc is the cyclotron frequency corresponding to the elec- 
tron pocket and l/r„ is the scattering rate of the elec- 
trons from the vortices in the mixed state. The analysis 
of Stephen also leads to a tiny shift of the relative fre- 
quency of the quantum oscillations in the mixed state, 
of the order of 10~^. Thus, there is enough leeway that 
even a very large error in this estimate will not affect our 
conclusions. 

The organization of the manuscript is as follows. In 
Sec. II we compute the ARPES spectral function and 
show that disorder can destroy the evidence of electron 
pockets. Section III is devoted to an exact transfer ma- 
trix computation of Shubnikov-de Haas oscillations and 
the effect of disorder on it. Section IV contains a discus- 
sion of scattering of quasiparticles of the putative normal 
state from the vortices in the mixed state. In Sec. V we 
briefly summarize the salient features of our work. 



II. SPECTRAL FUNCTION 

A. Hamiltonian 

The Hamiltonian of commensurate DDW order in 
terms of the fermion creation and destruction operators, 
cjj. and Ck, in the momentum space is 

Hi= (^kcj^Ck + ek+QC^+qCk+q) 

weRBZ 

+ ^ {iW]^clci,+Q + h.c), 

keRBZ 

where the ordering wave vector Q = (tt, tt), and ek is the 
single particle spectra. The lattice constant is set to be 
unity for simplicity. The reduced Brillouin zone (RBZ) 
is bomidcd by ky ± ~ ±7r. Wc define ek by 

ek = — 2t{cos kx + cos ky) + it' cos kx cos ky 

- 2t"{cos2kx+cos2ky) (2) 

and the DDW gap by 

Wo 

Ww = —{cos k^ - cos ky). (3) 
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B. Disorder 

Potential disorder in real space with a finite correlation 
length Id is modeled by 

nO^^/dxc'^GCx), (4) 

where the disorder averages are (G'(x)) — and 
(G(x)G'(y)) = (5(x — y); the disorder intensity is set by 
gy . This disorder Hamiltonian in the momentum space 
is then 

^2= 51 ^(ki,k2)4^Ck, (5) 

ki,k2e-BZ 

where the matrix elements are 

F(k,k + q) = |^e-^u(q), (6) 

and u(q) is 

w(q) = ^ /dy G(y)e-^'iy, (7) 

satisfying the conditions of (^(q)) = and (it(q)M(q')) = 
(5(q + q'). In practice, we generate u(q) directly with 
the desired statistical properties and then compute the 
matrix elements in Eq. ([5]). 

C. Computation of ARPES spectral function 

Once the full Hamiltonian H = Hi + H2 is generated, 
it is diagonalizcd by the transformation Ck = J2i ^k,i7i, 
where 7; is the annihilation operator of quasiparticles 
with energy Ei. The coefficients Pk.; and energy Ei are 
obtained through an exact numerical diagonalization pro- 
cedure. Finally, the ARPES spectral function A(k, uj) at 
a temperature T is given by: 

A{k,u;)^27rY,\P^i\^nid{Lu^Ei), (8) 

where n; = 1/[1 + exp{{Ei — fij/ksT)] is the fermion 
occupation number. Note that the numerical implemen- 
tation of Eq. ([5]) requires an approximation of the delta 
function by, for example, a Lorentzian distribution. 

We discretize the BZ with a mesh of size 80 x 80, and 
diagonalize the corresponding Hamiltonian. The param- 
eters we choose for YBCO at 10% doping are: t = 0.3 eV, 
t' = 0.3t, t" = t'/9.0, and Wq = 0.0825 eV, same as be- 
forcii The chemical potential 11 is set to be —0.2627 eV. 
These parameters yield a hole doping of Uh ~ 10%. The 
temperature T = 10 K is chosen, where a typical ARPES 
experiment is performed. For gv/i^Tv) = O.lt, the quasi- 
particlc life time for Id = is of the order r ~ 10~^^ s 
from Fermi's golden rule, which is a reasonable valuei^i^ 
Note that the band width is 8t. The scattering rate for 
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FIG. 1; The spectral function A(k,iL!) at a; = /i with cor- 
related disorder corresponding to Id = 4. The remaining 
parameters are stated in the text. 



correlated disorder for a finite Id ^ 4 is considerably 
smaller, as can be seen from Eq. IB] The final spectral 
function is obtained by averaging over 20 - 50 disorder 
configurations until no difference is detected upon fur- 
ther averaging. A typical result with disorder correlation 
length Id = 4, in units of the lattice constant, is plotted 
in Fig. [TJ The electron pockets are barely visible, resem- 
bling experimental observations. From Fermi's golden 
rule, the scattering rate is proportional to the square 
of the matrix element between the initial and the final 
states, y(k, k-l-q) ~ exp(— (7^/|,/2) (see Eq. (|6])), where 
q is the momentum transfer in the scattering process. On 
average, the scattering rate is therefore proportional to 
exp(— g^/|j/2); where qt is a typical momentum transfer. 
In particular, qt is roughly the size of a pocket. Because 
electron pockets are much smaller than hole pockets, the 
scattering rate is greater for electron pockets. The hole 
pockets centered at (±7r/2, ±7r/2) have vanishingly small 
spectral function on the back side due to the coherence 
factors and appear as Fermi arcs instead^^ whose lengths 
are further reduced by disorder. 

We also demonstrate the dependence of the spectral 
hmction ^(k, w = /i) on Id in Fig. [21 The disorder corre- 
lation lengths are Id = 0,2,8, 16 for panels (a) through 
(d), respectively. There are three distinct regimes de- 
pending on Id- For small Id, the electron and hole 
pockets will be almost equally scattered. As a conse- 
quence, the spectral function is smeared out everywhere; 
see Fig.[2l^a). For intermediate values of Id, for example 
= 2 in Fig. ^h) and Id = 4 in Fig. [1] scattering is 
more prominent for electron pockets, resulting in a pic- 
ture consisting of only four Fermi arcs. Finally, as the 
correlation length Id increases further, Fig. [2Kc) and (d), 
the electron pockets reappear. Indeed, though more dis- 
order scattering occurs on the electron pockets, the spa- 
tial variation of disorder, hence the net effect of disorder. 
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FIG. 2: The spectral function A(k, — fi) for gv/{2iT) — 
O.lt. The correlation lengths are Id = 0,2,8, 16 from panels 
(a) through (d), respectively. The remaining parameters are 
given in the text. 



becomes weaker. 
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FIG. 3: The energy dependence of the spectral function 
j4(k,aj) at two k points in the Brillouin zone, as indicated 
in the insert. 

To characterize the energy dependence of the spectral 
function, we compute A(k, to) as a function of u> at two k 
points in the Brillouin zone. One of them is at the inter- 
section TY line with the inner side of the Fermi surface, 
the other is situated on the electron pocket along the di- 
rection TAI; see Fig. [31 The disorder correlation length 
was chosen to he Id =4. Although there are peaks at 
uj fj, for both, the peak corresponding to the electron 
pocket is significantly suppressed by disorder; the sec- 
ond peak at w — /X ~ — 0.2eV is clearly an artifact of our 
simple theory and such high energy states would surely 
decay once correlation effects are taken into account by 
the creation of particle- hole pairs. 

In Fig.m A(k, oj) is plotted as a function of both k and 
u!. The horizontal axis is along the path T ^ Y ^ M ^ 
r in the Brillouin zone. The spectral function is negli- 



gibly small outside the reduced Brillouin zone bounded 
by kx ± ky ^ ztn^ and consequently there are no peaks 
in the central region of Fig. [H Close to k = (tt, 0) and 
Lu ~ 1^1, A(k, uj) has very small intensity due to long range 
correlated disorder, consistent with our previous observa- 
tion that the electron pockets are most likely unobserv- 
ablc. 
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FIG. 4; (Color online)The gray scale plot of the spectral func- 
tion ^(k, u)) as a function of both k and u). Red dashed curves 
indicate the quasi-particle dispersion relation. A horizontal 
dashed line shows the chemical potential /x. The remaining 
parameters are given in the text. 



III. SHUBNIKOV-DE HAAS OSCILLATIONS 

A. The transfer matrix method 

Let us now consider the effect of disorder on 
Shubnikov-de Haas (SdH) oscillations of the conductivity, 
axx- The tight-binding Hamiltonian on a square lattice 
in a sample of dimension N x M, with the lattice constant 
set to unity, is 



H 



' ■■CjCj -f- h.c. 



(9) 



where Ci is the fermionic annihilation operator at the site 
i. The spin degrees of freedom are omitted for simplicity. 
The hopping amplitude tij vanishes except for nearest 
and next nearest neighbors. To include two- fold commen- 
surate DDW order, the nearest neighbor hopping ampli- 
tudes are chosen to be 



(10) 



where (n, m) are a pair of integers labeling a site: i = 
nx -I- my, and Wo is the DDW gap; for the next near- 
est hopping tij = t' . The on-site impurity energy ei is 
defined by 



Ci = 



Z 



(11) 
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which is analogous to Eq. (|4]). To model correlated 
disorder, we set the disorder averages (Gr) = and 

{GrGr') = <5r,r', and Z = e^l''l^/^'o is a normalization 
factor. Vo parameterizes the disorder intensity. Note that 
Eq. lfTTj) reduces to a = Vod in the limit ^ 0, and ei 
becomes uncorrelated random variables. A constant per- 
pendicular magnetic field B is included via the Peierls 

phase factor aij = ^ // A • dl, where A = (0, —Bx, 0) 
is the vector potential in the Landau gauge. We note 
that a perpendicular magnetic field even as large as 60 T 
has little effect on DDW order 

In this section we choose t = 0.29 eV, t' = 0.1 eV, 
and Wq = 0.065 eV. The chemical potential is set to be 
^ = —0.28 cV. Note that these parameters are slightly 
different from those in the previous section, although the 
hole doping is again ^ 10%. We have left out the third 
nearest neighbor hopping, which greatly complicates the 
transfer matrix calculation without offering any particu- 
lar insight. The disorder intensity Vq = QAt leads to a 
quasi-particle life time of the order of ^ 10^^^ s in the 
limit of //J = 0. The magnetic field ranges from i? = 20 T 
to = 75 T, representative of the quantum oscillation 
experiments. The only relevant length scale here is the 
magnetic length Ib = y/K/eB, which for i? = 20 T is 
~ 15a, a being lattice constant equal to 3.85A. 

Now consider a quasi-lD system, N ^ M, with 
a periodic boundary condition along y-direction. Let 
^'n = {i^ns, ■0ri,2, • ■ • , ipnAi)^ bc tlic amplitudes on the 
slice n for an eigenstate with a given energy E, then the 
amplitudes on three successive slices satisfy the relation 
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(12) 

where iJ„ is the Hamiltonian within the slice n, and 
the matrix T„ corresponds to the hopping between the 
slices n and n + 1. Tn is tridiagonal, as electrons can 
hop from a site on slice n to three sites on the slice 
n -\- \. All postivc Lyapunov exponents of the transfer 
matrix^ 7i > 72 > • ■ • > 7m, are computed by iterating 
Eq. (|12p and performing orthonormalization regularly. 
The convergence of this algorithm is guaranteed by the 
well known Oseledec theorem.— For the above parame- 
ters a transverse dimension corresponding to M = 40 is 
sufficient. Equation was iterated 10^ to 10^ times 
until the relative errors of less than 1% of all the Lya- 
punov exponents were achieved. 

B. Computation of a^x 

The conductivity a^x at zero temperature is obtained 
from the Landauer formulai^^i^i^^ 

2 ^-f 1 

The SdH oscillations of Uxx at zero temperature for 
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FIG. 5: Fourier transform of SdH oscillations at zero temper- 
ature for /d = in arbitrary units. Insert shows oscillations 
as a function of the inverse magnetic field. Note the presence 
of higher harmonics. The parameters are described in the 
text. 



Zd = is shown in the insert of Fig.[5l A third order poly- 
nomial was used to subtract the background. The Fourier 
transform is shown in the main panel of Fig. [5] . Clearly, 
there are two main oscillation frequencies Fi — 490±30T 
and F2 = 915 ± 40T, corresponding to electron and hole 
pockets respectively, although a few harmonics are also 
visible. Though the first peak at Fi agrees with experi- 
mental observations, the second peak at F2 has not yet 
been observed. 

For correlated disorder, the Fourier transform of Uxx is 
plotted in Fig.[6l In all the cases, two main frequencies of 
Fi - 500T and F2 ^ 900T are prominent. As for A{k, uj), 
an increase of I d increases the amplitudes of Fi and F2 , 
because the effective disorder becomes weaker with the 
increased correlation length. Also the higher harmonics 
are more visible. There is, however, a sharp distinction 
between the dependence of the two physical quantities, 
which becomes clear when we consider the white noise 
case: A(k, oj) is completely smeared out in the momen- 
tum space because white noise scatters between all pos- 
sible wave vectors; sec Fig. [3] (a). The coherence factors 
are of crucial importance for the spectral function. In 
contrast, the SdH oscillations are damped by the Dingle 
factor, which is parametrized by a single lifetime and dis- 
order enters in an averaged sense. This striking contrast 
is clear in Fig. [5] The surprise is that impurity scatter- 
ing affects the electron pocket more than the hole pocket, 
which is remarkably robust even for white noise. There 
must be a piece of physics missing, if we are to explain 
why the hole pocket is not observed in SdH. This miss- 
ing physics we argue is the vortex scattering rate which 
affects the two pockets very differently. 



IV. VORTEX SCATTERING RATE IN THE 
MIXED STATE 

We have shown previously^ that the scattering rate of 
the DDW quasiparticles, corresponding to the electron 
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FIG. 6: Fourier transform of SdH oscillations for corre- 
lated disorder at zero temperature. The correlation length 
Id = 1, 2, 4, 8, 16 for panels (a) through (e), respectively. The 
remaining parameters are described in the text. 



pockets, from vortices in the mixed state is given by. 

1 \ _ A2 / TT 



(14) 



where the cyclotron frequency lOc is determined from the 
band mass mi, to be (restoring the lattice spacing a here) 



2mb 



2t' + U" 



4 



(15) 



Here A^, a measure of the amplitude of the supercon- 
ducting gap, is an average over the disordered configu- 
ration of the vortices, which is likely to be insensitive to 
the symmetry of the order parameter. Similarly, we have 
shown that the vortex scattering rate for the nodal DDW 
quasiparticles at the hole pocket is given by (noting a ty- 
pographical error in Ref. 
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A2 



27r 



(16) 



The relevance of nodal quasi particle dispersion is also 
pointed out in Ref. ? We have introduced a character- 
istic mass and a frequency scale by (nodal quasiparticle 



are actually massless) 



m 



VpVD 

eB 



(17) 
(18) 



The reason for this is that these are precisely the scales 
that enter a SdH calculation of the nodal quasiparticlesi^ 
It can be seen from the above equations that the vortex 
scattering rate for the nodal particles is proportional to 
I/^/Wq. This idi verges as Wq — s- 0. But this is not 
unphysical, as, in that limit, there is a phase transition, 
where the Fermi surface reconnectsJ^ The nodal fermions 
are obviously no longer a valid description. However, our 
mean field theory cannot be trusted to predict the precise 
scattering rate at this quantum critical point. We have 
restricted ourselves to a regime that we believe is deep 
inside the broken symmetry phase where our treatment 
should be a good guide, and where we believe that the 
present experiments are performed. A proper treatment 
of this scattering rate at the point where Wq collapses is 
an open question that may shed light to the physics of 
the striking strange metal phase. 

The physical picture underlying the calculation of 
Stephen^S is that the lifetime obtained from the imagi- 
nary part of the self-energy corresponds to a situation as 
if the vortices are static impurities. Here vp ^ 2^/2at/h 
is the magnitude of the velocity in a direction normal to 
the hole pocket and vd = Woa/-\/2S is the velocity in a 
direction orthogonal to it; neither t' nor t" enter in the 
leading order. Therefore, 
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3/2 



1 



(19) 



The cancellation of is interesting but cannot of course 
be an exact result considering the approximations in- 
volved in Stephen's analysisi^ The correctness of which 
relies entirely on the formulation of the averaged Green 
function of Stephen, where this average is carried out 
in the real space, and the nature of the normal state 
(whether or not it contains electron and hole pockets in 
the momentum space) does not enter, as the normal state 
Green function does not contain the superconducting or- 
der parameter. On the other hand, we had roughly esti- 
mated JTifc « 1.27me and m* « 2.72TOe for a typical set of 
parameters)^ where me is the free electron mass. Thus, 
it is reasonable that 
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1 



UJcTv 



(20) 



which should be robust with respect small changes of the 
band parameters and the DDW gap. This would lead 
to a strong suppression of the oscillations corresponding 
to the hole pocket as compared to the electron pocket, 
as the Dingle factors, I'e,/! Siie exponentially sensitive to 
the product of the cyclotron frequency and the vortex 
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scattering lifetime: Ve = e-"" l^'^--^-')- and Vh = Vf^ . 
The previously estimated lifetime of the electrons;^ 

(^■^^ = 1.5 X lO^^s-i, (21) 

for B = 40 T and 3^2 = 60 T should be a rough guide. 

The whole analysis is predicated on the assumption 
that the quantum oscillation frequencies are unshifted 
from the putative normal state (the DDW state in this 
case), for which we now provide some support from the 
analysis of Stcpheni^ see, however Ref. |3^. The for- 
mula for the relative frequency shift in terms of physical 
parameters (absolutely essential because they are all ef- 
fective parameters) is 

^^1^ (22) 

where A(i7) is the zero temperature superconducting gap 
in a magnetic field; for a d-wave superconductor we have 
to use an appropriate average, as above. The factor hiOc 
for free electrons is given by 

hujc = 1.34 X IQ-^ X H{Tes\a)eV. (23) 

As will be seen below, it will make little difference if 
we use the mass determined from experiments. The 
magnetic field ranges between 30T and 65T. Let's take 
H = 40 T. Then 

ftwc = 5.36 X IQ-^eV. (24) 

We know little about the zero temperature gap, especially 
in the underdoped regime, where there are fluctuation 
effects. As the simplest assumption, we use BCS mean 
field theory: 

A{H) ^ A{0)^1- H/Hc2. (25) 

and 2A(0) = 3.52 Tc, which results in 

A(0) = 8.6 meV (26) 

for Tc = 57.5 K.i^ For simphcity we choose A(0) w 
10 meV, and as a rough guide from experiments,—— Hc2 ~ 
60 T, although it will make little difference even if it 
were 100 T. For 10% doping we need a chemical potential 
/i w -0.26 eV. Thus, we get 

AF 

— = 4.6 X 10"^. (27) 

Even if this estimate were incorrect by several orders of 
magnitude, our conclusions that the frequency shift is 
negligible should be valid. 

Consider the white noise as an example; see Fig. [5] We 
first filter out the two peaks, invert the Fourier transform, 
and then multiply by the Dingle factors corresponding to 
the vortex scattering rates discussed in this section. Us- 
ing Eg . [20l and Eg. [21] and Fourier transforming back this 
procedure essentially wipes out the peak corresponding 
to the frequency of hole pocket in the resulting transform, 
as shown in Fig. [T) 
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FIG. 7: Fourier transform of SdH oscillations shown in Fig. [5] 
after taking into account vortex scattering rates in the mixed 
state, as discussed in the text. Note that the amplitude of 
the peak corresponding to the hole pocket at about 900 T is 
essentially wiped out. 



V. CONCLUSIONS 

In summary, we have shown that two of the most puz- 
zling features in underdoped cuprates, the nonexistence 
of electron pockets in ARPES and the lack evidence of 
the hole pockets in quantum oscillations can be plausi- 
bly resolved by considering correlated disorder and the 
vortex scattering rates in the mixed state. With respect 
to the latter, it is crucial that the nodal fermions are 
described by a relativistic spectrum while the quasipar- 
ticles corresponding to the electron pocket are described 
by a non-relativistic spectrum determined by the bottom 
of the band. 

It is important to determine if the calculation of 
Stephen can be so closely taken over after modifying for 
the presence of nodal fermions. We believe that it can be 
because the relative frequency shift turns out to be very 
small, of the order of 10~^. Thus, even if the estimate 
is off by some orders of magnitude, we still have enough 
leeway. This argument has been recently challenged^ in 
a calculation of the density of states, however. Clearly, 
further investigations will be extremely valuable and are 
in progress. 

One of the surprising conclusions of the present work, 
which we believe is on firm grounds, is that even very 
strong disorder, such as white noise, has only a modest 
effect on quantum oscillations, while it has a much larger 
effect on the ARPES spectra. This is because ARPES 
spectral function depends on the coherence factors, which 
act as Wannier functions that are naturally very sensitive 
to disorder. In quantum oscillations disorder appears 
primarily as an averaged lifetime in the Dingle factor. 
Elastic scattering has little direct effect on the Onsager 
quantization condition. 

It is also a firm conclusion of our work that elastic 
disorder scattering from impurities cannot be responsi- 
ble for wiping out the hole pocket frequency, while keep- 
ing the electron pocket frequency more or less intact. 
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In fact, quite the opposite seems to be true. Thus, 
the vortex physics in the mixed state appears to be of 
paramount importance, especially the contrast between 
the non-relativistic electrons around the anti-nodal point 
and the relativistic nodal fermions at the nodal points of 
the Brillouin zone. 

We also note an interesting paper— regarding an anal- 
ysis of various masses involved, which take into account 
electron-electron interactions. In the future, it would be 
interesting to pursue such an analysis of residual electron- 



electron interactions in the present context. 
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